Interaction-driven effects on two-component Bose-Einstein condensates 
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We investigate the role of the interparticle-interaction strength in the distribution of two species 
of atoms inside a condensate. We focus upon the study of systems for which the minima of the 
trapping potentials for the species are displaced from each other by a distance that is small compared 
to the size of the total condensate. We show that in a small range of the interparticle-interaction 
strength the distribution of species undergoes dramatic changes, and exhibits a variety of different 
features. We demonstrate that this behavior can be easily understood in terms of the Thomas-Fermi 
approximation. This effect may be useful in experimentally determining the values of the scattering 
lengths. 
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Much experimental and theoretical work in Bose- 
Einstein Condensation deals with systems composed of a 
mixture of two distinct species of atoms [[j} |[ [| [| |[ || . 
For example, one mixture commonly used is that of 
atoms of 87 Rb in two different hyperfine states \F — 
l,m/ = —1) and \F = 2,m/ = 1) This mixture has 
shown itself to be very useful for experimentally gener- 
ate topological modes, such as different types of vortices 
j^, |J, |5| ■ Experimentally these mixtures are under condi- 
tions such that the species behave as "effectively distin- 
guishable" , and have been observed to separate partially 
in space An interesting issue to investigate is how 
the particles distribute themselves inside the condensate 
depending on the interparticle interaction, specially be- 
cause there still exists uncertainty about its numerical 
value. 

From a theoretical point of view, H. Pu and N. P. 
Bigelow studied the ground state properties of differ- 
ent mixtures in a spherically symmetric trapping poten- 
tial, fn Ref. 0, they displayed the density profiles of 
both components for various combinations of scattering 
lengths below phase separation, testing the validity of the 
Thomas- Fermi approximation. Later, E. Timmermans 
extended the study to phase separation || . In particular 
he showed that in a spherically symmetric trap the less 
repulsive component remains inside a sphere while the 
other component lies in a spherical shell around it. 

In this work we analyze the case when one of the 
species is in a slightly shifted potential with respect to 
the other species, and thus the spherical symmetry is 
broken. We find that by varying the interparticle in- 
teraction strength by less than 10%, the particles rear- 
range themselves in very different configurations inside 
the condensate. As far as we know, this problem was 
only studied for a particular set of scattering lengths |j] 
and it still lacks an easy interpretation. Here we show 
that the distribution of species can be easily understood 
in terms of the Thomas-Fermi (TF) approximation. We 
compute the exact Gross-Pitaevskii (GP) solutions for a 
large number of particles and find that TF predictions 



describe with great accuracy the geometrical properties 
of the distribution of species within the condensate. 

In order to describe the wave functions of two-species 
condensates one has to solve the coupled Gross-Pitaevskii 
equations Jl0| : 
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where Ni and Mi denote the number of atoms and the 
mass of the species i, respectively. The factors Gk,i — 
Uk,iU, with Uk,i being the relative interaction strengths 
between species k and I, and Vi the potential seen by 
species i. Hereafter, we consider a two component system 
and set the most repulsive component |1) in the "Fi state, 
fixing ui.i = 1 (> 1*2,2), so U — Ann a/Mi, a being 
the scattering length of the |l)-species. In addition, for 
simplicity, we set Mi = M2 = M. 

Qualitative analysis — In order to study qualitatively 
how the species rearrange themselves inside the conden- 
sate when varying the relative interparticle strength 1*1,2, 
we make use of the Thomas Fermi (TF) approximation, 
which consists of neglecting the kinetic energy, and thus 
removing the laplacian terms in Eqs. We will con- 
sider a system in which both components are in an axially 
symmetric trap and the |2) component has the minimum 
shifted in the z direction by a v alue — d. We make a 
change of variables according to \J (M/2)uj z f ^ r, where 
ui r and u z are the trap angular frequencies in the r and 
z coordinate, respectively. Futhermore, by defining the 
aspect ratio A = u> r /u> z , the potentials Vi, written in 
cylindrical variables, read 



AV + z 2 



and 



v 2 = av + { z + dy 
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The sign of the determinant A = 1*1,11/2,2 -«i 2 de- 
fines two different features for the distribution of par- 
ticles. When A > and 1*1,2 is small, there exists a 
large coexistence region. As 1*1,2 is increased this region 
decreases. At A = 0, a phase separation takes place re- 
ducing the coexistence region to an interface. We shall 
consider the cases A > and A < separately. 

Coexistence A > — The solution of Eqs. ([!]) in 
the TF approximation can be easily obtained and has 
the following different expressions depending on whether 
there exists any overlap between the wave functions of 
both species. 

a) In the region where only one wave function is non 
vanishing ( |*,| 2 ^ and |* fc | 2 = 0, for i ^ k) the TF 
equations are decoupled and the solution reads 
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|* 2 r = ^2-AV-(z + d) 2 ]/(G 2 , 2 7V 2 ). (4) 

b) In the region where both wave functions are non 
vanishing l^l 2 > 0, the solution may be written after 
some algebra as: 
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ui. 2/1^2,2, 02 = "1,2/1/1,1 and Bi 



where 0\ 
0i)/{UNi&). 

The surfaces 1S1 and 5 2 , defined by equating the ex- 
pression inside the square bracket in |*i| 2 and |* 2 | 2 of 
Eq. (^J) respectively to zero, determine the boundary of 
the coexistence region. It may be seen that these sur- 
faces are ellipsoids centered in d\ — 0id/(l — 0i) and 
d 2 = — d/(l — 02) along the z-axis, respectively. It is 
worth mentioning that this result does not depend on 
any other quantity, as for example the number of par- 
ticles or the frequencies of the trapping potential. The 
factor 1 — 0i in the denominator makes these displace- 
ments diverge when 0j is close to unity, and one can 
guess that for these values, some dramatic effects in the 
redistribution of particles could take place, as we shall 
discuss later. 

Phase separation A < For these interac- 
tion strengths, in the Thomas-Fermi approximation, the 
species are confined to two separated regions within the 
condensate. The boundary surface S s between these two 
regions may be obtained by equating the pressure (TfJ 



on both sides of the interface, yielding 



(7) 
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Assuming that the wave functions on each side are given 
by the expression ([!]) and defining a = \J G 2j2 /Gi,i, the 
interface obtained is 
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with 
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This surface is an ellipsoid whose center is displaced in 
d s = —d/ (1— a) along the z-axis. Note that the quantities 
R s and d s do not depends on iti j2 . It is easy to prove that 
for A = the surfaces verify Si = 5 2 = S s . 

In what follows we consider only A = 1 in which case 
the surfaces Si are spheres with radii Ri. 

In the TF approximation, the density contours pt = 
|*j| 2 (i — 1, 2) inside the coexistence region are spherical 
surfaces with radii R < R4 ( R > R4) if 0i < 1 (0i > 1) 
centered in dj. 

Numerical results — On the one hand, we compute 
the displacements di and radii Ri in the Thomas-Fermi 
approximation. Note that these radii depend on the 
chemical potentials, and thus also on both the number of 
particles and the trap frequencies. On the other hand, in 
order to obtain the exact densities, we solved the Gross- 
Pitaevskii equations, using a steepest descent method. 

In particular, we used the relative intraparticle inter- 
action strength of atoms of 87 Rb in the two different hy- 
perfine states given, within a 1.2% of error, in Ref. [Q, 
141,1 = 1 and it 2j2 = 0.94. We chose a trapping potential 
with an angular frequency oj r — 2tt x 7.8 Hz, and set the 
number of particles of each species N% = iV 2 = lxlO 7 . 

In Fig. |l] we show the Gross-Pitaevskii density con- 
tours pi in the y = plane, together with the Thomas- 
Fermi displacements di and radii Ri, for d — 0.4 /zm. 
From the first to the last row we vary 1/1,2 by less than 
10%. For i/i )2 = 0.93, which corresponds to the first row 
of the graph, it can be seen that a large coexistence region 
still exists. This region is given by the intersection of the 
circles defined by Ri and di, which leaves outside only a 
little region on the top of the condensate filled with |1) 
type particles. It can also be seen that inside the coexis- 
tence region, the GP density contours are well-described 
by spherical surfaces centered in the points d\ and <i 2 for 
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FIG. 1: GP density contours of pi and p2 are shown in the 
first and second column, respectively. From top to bottom of 
the graph, the rows correspond to relative interparticle inter- 
action strengths ui, 2 = 0.93, 0.94, 0.96, V0.94 and f. The 
dashed, dotted and dash-dotted lines correspond to TF radii 
Hi, i?2, and R s , respectively. The dot, cross and triangle 
indicate di, d,2, and d s , respectively. The trapping potential 
displacement is d = 0.4/im. 



P\ and P2 respectively, as predicted with the TF approxi- 
mation. For the second row we consider 1*1,2 = 0.94, and 
with this value d\ and R\ diverge. The density contours 
inside the coexistence region are quite planar surfaces 
for the |1) component, as expected from the TF anal- 
ysis because of the above-mentioned divergences. For 
1*1,2 = 0.96, which is displayed in the third row, three 
phases already exist: pure |1) and |2) components and a 
coexistence region. Once more, the contours seem to be 
in agreement with formulae (|j) and (|^) with only a small 
departure in the region next to the boundaries. 




FIG. 2: Thomas- Fermi displacements (thick dashed 
line), lefel (thin dashed line), radii R\ (thick solid line) and 
R2 (thin solid line), and the TF radius of the condensate Rtf 
(dash-dotted line) as function of the interparticle interaction 
strength ui,2 for the same conditions as in Fig. 1. 

Phase separation occurs for u\$ — \/0.94, which corre- 
sponds to the fourth row, and thus the coexistence region 
is reduced to a surface. It is evident that although the 
GP solutions exhibit some overlap the general feature of 
the condensate is well-described by the TF approxima- 
tion. Of course, if we take a smaller number of particles 
into account, the overlap will be greater. The last row 
corresponds to 1*1.2 = 1 which is quite similar to the pre- 
vious figure. This is due to the fact that formulae (||) 
and ([iol ) do not depend on ui, 2 , as we stated before. 

In Fig. H we display the TF quantities, d, and R4 
together with the radius of the condensate -Rtf as func- 
tions of 1*1,2- For 1*1.2 < 0.9, all these quantities are 
smooth and monotonous functions, while in the inter- 
val displayed in the figure, the quantities related to |1) 
species exhibit an abrupt behavior. This effect suggests 
that if one wants to test experimentally the validity of the 
values of the interaction strengths within this interval by 
determining the displacements di, the error in their de- 
termination should not affect the desired quantities too 
much. 

As a final illustration, in Fig. || we display the density 
contour for two different sets of relative scattering lengths 
to describe the interaction of the same two species of 
87 Rb. One set, which is assumed to be more accurate 
§, is 1*1,1 = l,«i,2 = 0.97 and 1*2,2 = 0.94. However, to 
tell the truth, up to our knowledge, the error of 1*1,2 is 
not given anywhere. For these values of the interaction 
strengths the determinant is neg ative (A = -9xl0" 4 ) 
and thus the system is phase separated, and hence the 
coexistence region in the Thomas Fermi approximation 
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FIG. 3: The first and second row show the GP density con- 
tours of pi and p2, respectively. From the first to the third 
column we consider d = 0.4, 1, and 1.5 /im. The dashed lines 
correspond to the set = l,ui,2 = 0.97 and ti2,2 = 0.94 
while the solid lines correspond to the relative interaction 
strengths ui.i = l,m,a = 108/109 and u 2 , 2 = 108.8/109. 
The TF points d 3 (for the first set of parameters), and di and 
di are marked with a triangle, a dot and a cross, respectively. 

is reduced to an interface. The other, previously used, 
set of parameters H is u\ t \ = 1, 1*1,2 = 108/109 and 
1*2,2 = 108.8/109, and verifies A =' 1.6xl0~ 2 . Since 
the determinant is positive a coexistence region exists. 
Moreover, j3\ — 0.99 is less than 1, and the contours 
of the |1) component at the coexistence region are still 
concave upwards. 

We used three different displacements d = 0.4, 1, 1.5 
/xm for the trapping potentials. For the first set of inter- 
action strengths the corresponding TF displacements are 
d s = —13, —32, —49 /im, respectively. While for the sec- 
ond set, the displacements are d\ = 54, 135, 202 (im and 
c?2 = —44, —109, —163 /im. It may be seen that the dis- 
tribution of particles is very different in both cases, espe- 



cially for small displacements, and this fact can be used 
for testing experimentally the values of the scattering 
lengths. It is worth mentioning that in the experimental 
conditions of Ref. ||, as they obtain integrated densi- 
ties, one cannot distinguish between a phase separated 
system (A < 0) and a system in which the components 
are overlapped (A > 0). 

Finally, having in mind the relations j3\ = d\/ {d + d\) 
and (32 = {d + g?2)/<^2, we wondered wether having the 
experimental density contours in a given plane, say y = 
0, one can accurately determine the scattering lengths 
(see definitions of Pi) by estimating the displacements 
d^ In order to answer this question we have done the 
following test. We have used the information of the GP 
density contours assuming that they represent the exact 
experimental data. For a given density contour z{x) in 
the y = plane we computed the slope b = dz/dx at 
each point (x, z) of the curve. Then we calculated the 
intersection point between the line perpendicular to the 
density contour and x = 0, which gives the center of the 
spherical surface. This point is z c = z + x/b. Plotting 
z c for points all over the |1) (|2)) component contour we 
found well-defined plateaus at z c — ( z c — ~ d) for 
points outside the coexistence region, and in z c = d\ ( 
z c = d^), where d* are the estimates of the TF d t . We 
used this information to determine the scattering lengths 
and found an error in their determination of about 0.1%. 
Note that, even considering an error in z c of about 30 
/tm and using the expressions of the uncertainty in 
the value of 1*1,2 turns out to be below 0.5%. A similar 
procedure could be carried out with experimental data. 

In summary, we show that the distribution of com- 
ponents is strongly ruled by the interparticle interaction 
strength, and for the number of particles we have consid- 
ered, the way the particles rearrange inside the conden- 
sate can be easily understood in terms of the Thomas- 
Fermi approximation. We also outline a possible proce- 
dure to experimentally determine the relative scattering 
lengths. 
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